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We describe one-dimensional photonic crystals that support a guided mode suitable for atom 
trapping within a unit cell, as well as a second probe mode with strong atom-photon interactions. 
A new hybrid trap is analyzed that combines optical and Casimir-Polder forces to form stable traps 
for neutral atoms in dielectric nanostructures. By suitable design of the band structure, the atomic 
spontaneous emission rate into the probe mode can exceed the rate into all other modes by more 
than tenfold. The unprecedented single-atom reflectivity ro > 0.9 for the guided probe field should 
enable diverse investigations of photon-mediated interactions for ID atom chains and cavity QED. 
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New opportunities in Atomic, Molecular, and Optical 
Physics and Quantum Information Science emerge from 
the capability to achieve strong radiative interactions 
between single atoms and the fields of nanoscopic op- 
tical waveguides and resonators [1 . For example, strong 
atom-photon interactions in lithographic structures [2]-^ 
could be used to create quantum optical circuits with 
long-range atom-atom interactions mediated by single 
photons [71 H]. Moreover, linear arrays of atoms ra- 
diatively coupled to nanophotonic waveguides exhibit a 
wide range of remarkable phenomena, including coher- 
ent transport of atomic emission [9 11 , guided superra- 
diance and polaritons p!2HT4] , as well as highly reflecting 
atomic mirrors [TF, T6^ . The interplay of atomic emission 
into the waveguide and photon-mediated forces can lead 
to self-organization of atoms into exotic spatial configu- 
rations along the waveguide J7l [18] . 

A long-standing obstacle to this scientific frontier is 
the challenge of trapping atoms in vacuum near dielec- 
tric surfaces ('^ lOOnm) while at the same time achieving 
strong interactions between one atom and photon. A far- 
off resonance dipole-force trap (FORT) [19] can provide 
atomic localization by using modes of the dielectric for 
optical trapping [ 2QH22] and has been used to trap cold 
atoms within hollow-core optical fibers [23H26] and ex- 
ternal to fiber-taper waveguides [27H29] . 

Motivated by these advances, in this manuscript we 
present principles for the design of optical traps and 
strong atom-photon interactions in one-dimensional (ID) 
photonic crystal waveguides. We analyze the potential 
[/tot(r) due to light-shifts from a FORT [3QH32] together 
with Casimir-Polder (CP) interactions with the dielec- 
tric [33H36] (Figs. [TJa, e)). Despite the proximity of the 
surfaces, stable potentials /7tot(r) are achieved for mod- 
est optical intensities (~ 5mW//im^) for blue-detuned 
FORTs operated at a 'magic' wavelength for the D2 line 
of atomic Cesium [19]. A new possibility for trapping is 
also identified for which vacuum forces from CP interac- 



tions are exploited to close the trap perpendicular to the 
plane of structure, which would otherwise be unstable 
with either the FORT or the CP potential alone [37 . 

In addition to the waveguide trapping proper- 
ties, strong near-resonance atom-photon interactions of 
trapped atoms are found to arise for waveguides with 
properly tailored band structure [38ti42] . For practically 
realizable structures, we find ^id h' ^ 10, where ^id 
is the atomic decay rate into the (guided) probe mode 
and 7^ the rate into all other modes. One atom trapped 
within the structure could thereby attenuate a resonant 
probe with transmission |1 — rop ^ 10~^ [9, 16 . 

As illustrated in Fig. [l] we focus on two of the simplest 
quasi- ID photonic crystal geometries. The first waveg- 
uide consists of a single silicon-nitride nanobeam (refrac- 
tive index n = 2) with a ID array of filleted rectangular 
holes along the propagation direction; atoms are trapped 
in the centers of the holes (Fig. [ija)). The second waveg- 
uide consists of two parallel silicon nitride nanobeams, 
each with a periodic array of circular holes, with atoms 
trapped in the gap between the beams (Fig. [ije)). 

The design of a lI)-photonic crystal waveguide with 
distinct modes for optical trapping and strong atom- 
photon interactions is constrained by the region of the 
optical band structure containing a continuum of un- 
guided optical modes (i.e., the light cone indicated in 
gray in Fig. [TJb, f)). Modes within the light cone can 
still have large amplitude in the structure but radiate en- 
ergy into the surrounding vacuum leading to unaccept- 
able loss. The top of the vacuum light line is at the 
Brillouin zone boundary (X-point, where kxd = tt), so 
the lattice constant a is constrained by a < A/2, where A 
is the smaller of the (vacuum) wavelengths for trapping 
and probe fields. Here, kx is the Bloch wavevector along 
the waveguide axis x. 

Once a is fixed, additional guided modes can be 'pulled' 
below the light line by increasing the width and thickness 
of the structure. With appropriate modes below the light 
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FIG. 1: a) Schematic for the single nanobeam structure with 
dimensions {a,w,t, ha:, hy) = (367, 845, 825, 246, 745)nm. b) 
Band diagram for the single nanobeam in a) showing only 
bands with even vector symmetry about the y and z sym- 
metry planes. The trapping and probing bands are shown as 
thicker lines, with the trap ujt 1"^^ (probe uja/'^t^) frequency 
as a blue (red) dashed line, c) Field intensity of the blue 
trapping mode and d) field amplitude of the probe mode in 
the center plane z — ^ for the single nanobeam in a). Green 
spheres mark the locations of minima of the trapping poten- 
tial, e) Schematic for the double nanobeam structure with 
{a,w,t,d,g) = (335,335,200, 116, 250)nm. f) Band diagram 
for the double nanobeam in e) displaying only modes of even 
vector symmetry in z. The proximity of the two nanobeams 
results in a band structure composed of even (green) and odd 
(magenta) superpositions of single nanobeam modes. We fo- 
cus on the even parity supermodes due to their large field 
amplitude in the gap. g) Field intensity of the blue trapping 
mode and h) field amplitude of the probe mode in the center 
plane z = of the double nanobeam in e). The black dia- 
monds in b) [f)] mark resonances for finite structures of 81 
unit cells from Figs. [4] [s]. 

line for probing and trapping, the spacing of these modes 
at the X-point can be tuned by altering the size of the 
holes, which enables the probe mode to be resonant with 
the frequency uja of the atomic transition while simul- 
taneously matching the optical frequency ujt of the trap 
mode to a 'magic' frequency for the atom [T9j . 

Within this general context, here we consider only 
blue-detuned FORTs for which the trapping mode has an 
intensity minimum at the trapping site [43^. Our analy- 
sis is for the D2 line of atomic Cs with probe wavelength 
near the atomic resonance = 852nm and with a blue- 
detuned FORT at the magic wavelength At = 793nm 
[l9] . Note that our results are readily transcribed to 



other atomic transitions by way of the scale invariance of 
Maxwell's equations [44 . 

The photonic crystals are assumed to be suspended in 
vacuum and composed of SiN. Band structures are cal- 
culated using the MIT photonic bands software package 
j45] . Field profiles for guided modes are calculated using 
the finite-element-method (FEM) simulations [46]. Re- 
sults for the single and double nanobeam structures are 
presented in Fig. [l] 

With suitable guided modes for trapping in hand, we 
have developed numerical tools for evaluating the FORT 
and CP potentials inside the waveguide, and hence the 
total potential [/tot = ^fort + ^cp- The adiabatic po- 
tential [/FORT(r) is readily calculated using the electric 
field distribution of the trap mode, E(r) = u/e^(r)e*^^^ 
j3QH32] . Here Uk^ (r) is the periodic Bloch wave function 
at propagation constant kx. 

The surface potential Ucpi'^a) is determined from the 
formalism in Ref. for the imaginary component 

of the scattering Green's tensor Gsc(ra, r^, cj), which is 
the Green's tensor from Maxwell's equations for a point 
dipole at the atomic location with the vacuum con- 
tribution (i.e., no dielectric structure) subtracted. We 
evaluate Gsc(ra, r^, cj) numerically by adapting the pro- 
cedures from Ref. [36 , as described in [47 . 

Figures [2] and [s] display numerical results for /7cp(r), 
^FORT(r), and [/tot(r) for the single and double 
nanobeams for kx below the X-point. The calculations 
are for the 65'i/2, F = 4 hyperfine ground state of Cs for 
the FORT modes indicated in Fig. [ijc, g) [48]. For these 
initial calculations, we make the reasonable assumption 
for SiN that the dielectric constant e is frequency inde- 
pendent, e(r,co') e(r). 

/7tot(r) for the single nanobeam in Fig. |2] reveals that 
modest optical intensity is sufficient to overcome the at- 
tractive CP interactions and create a stable potential 
minimum in the center of the vacuum space at Vmin = 
within a unit cell. An atom would be localized at dis- 
tances (dx^dy) = (123,373)nm from the walls of the di- 
electric. The trap oscillation frequencies for a Cs atom 
would be {fx, fy, fz) ^ (612, 180, 484)kHz. 

For the double nanobeams [49], the FORT alone is 
insufficient to trap the atom, as the mode has a (weak) 
local intensity maximum along the z direction that repels 
an atom. However, the CP potential /7cp(r) along z 
provides the force necessary to overcome the repulsive 
optical force and to form a stable trap. The result is a 
hybrid optical- vacuum trap that circumvents the 'no-go' 
theorem for vacuum trapping alone [37 . 

Potentials for our hybrid trap are illustrated in Fig. 
[3] At the trap minimum Vmin = 0, an atom would be 
localized at distance dy = 125nm from adjacent surfaces 
of the dielectric beams. Oscillation frequencies for a Cs 
atom would be (/^, /^) (1013, 390, 57)kHz. 

As concerns strong radiative interactions, our struc- 
tures trap an atom in a region of large amplitude for the 
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FIG. 2: Trapping potentials for the single nanobeam structure 
in Fig. ^a) for Cs 65'i/2,F = 4 level and At = 793nm. 
(a) Casimir- Polder potential [/cp(r) and (b) total potential 
t/tot(r) = t/cp(r) + t/FORT(r) in the central z = plane, (c- 
e) show line cuts of Ucp (red solid), [/fort (blue dashed), 
and [/tot (blue solid) along the (c) x-, (d) y-, and (e) z-axis. 
Average trap intensity for a unit cell is 4.9mW///m^. 

probe field, leading to small mode volume per unit cell 
[47]. It is well known that atom-photon interactions can 
be further enhanced near a band edge j38H42] , where the 
density of states diverges due to a van Hove singular- 
ity. To quantify the radiative coupling, we determine the 
decay rate 7tot for a point dipole located at = for 
a structure with N unit cells [50^. FDTD calculations 
are performed to evaluate the classical Green's tensor 
G(ra, Ta^oj) and thence 7tot following Refs. [34l |35j [51], 
as described in [47 . 

Figures 4(a) and 5(a) display the diagonal components 
Im[Gii{i'd)] of the Green's tensor as functions of dipole 
frequency Vd = uJd /27r and relate to the emission rate of 
resonant point dipoles polarized along the i = or 
2;-axis for the single and double nanobeams. Firstly, in 
Fig. 4(a), Im[G xx] is enhanced along the (periodic) di- 
rection across a broad frequency range, and is suppressed 
in the y- and z- directions, as can be explained by the ori- 
entation of the induced array of image dipoles along the 
single nanobeam. When the source dipole is polarized 
along the x-axis, the image dipoles line up head-to-tail 
and, just below the X-point, constructively interfere to 
enhance dipole emission; y^ z orientations render destruc- 
tive interference and suppressed emission. Comparable 
suppression is not apparent for the double nanobeams in 
Fig. 5(a) since no such array of image dipoles is formed. 

Secondly, Figs. 4(a) and 5(a) display a series of res- 
onant peaks due to strong emission into various guided 
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FIG. 3: Trapping potentials for the double nanobeam struc- 
ture in Fig. |l|;e) for Cs 65'i/2, F = 4 level and At = 793nm. 
(a) Casimir- Folder potential t/cp(r) and (b) total potential 
t/tot(r) = t/cp(r) + t/FORT(r) in the transverse x = plane, 
(c-e) show line cuts of Ucp (red solid), [/fort (blue dashed), 
and [/tot (blue solid) along the (c) x-, (d) y-, and (e) ^-axis. 
Average trap intensity for a unit cell is 3.5mW//im^. 

modes. In the region near the Cs D2 line (i.e., Ud — i^a = 
ua/^tt = 352 THz), we find peaks in Im[Gxx] for the sin- 
gle nanobeam and in Im[Gyy] for the double nanobeam. 
These peaks are due to emission into our designated 
probe modes for the respective structures, where for the 
single (double) nanobeam(s), the probe mode is princi- 
pally polarized along the x-{y-) axis. Each peak is from 
a discrete set of propagation constants ki^"^ ~ irn/aN 
imposed by the boundary conditions for the finite struc- 
tures. Here, N is the total number of cells in the single 
(double) beam, and n < TV is an even (odd) integer. We 
find excellent agreement between the frequencies of these 
resonances and the band diagram of the probe mode ('di- 
amonds' in Figs. 1(b) and (f)) for various values of (n, N). 

The peaks become larger and narrower as /ci^^ ap- 
proaches the X-point, owing to the diminishing group 
velocity [40l |42l |47]. Beyond the X-point, the probe 
resonances disappear, leaving a broad background cor- 
responding to coupling into lossy (radiation) modes. 

On an expanded frequency scale around Ud — ^A^ 
Figs. 4(b) and 5(b) show calculated atomic decay rates 
7tot for the QP2,/2-,F' = 5 ^ ^Sii2-,F = 4 transition 
in atomic Cs [47 . When the atomic dipole is aligned 
along the principal polarization of the designated probe 
mode (x for the single beam and y for the double beam) , 
the emission rate ^id into the probe mode is strongly 
enhanced at frequencies corresponding to near the 
X-point. Specifically, for Vd = ^a large enhancements 
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FIG. 4: Green's tensor and total atomic decay rate 7tot versus 
source dipole frequency Ud for the single nanobeam at r = 0. 
(a) shows the diagonal components of the Green's tensor, 
Im[Gxx] (solid black), Im[Gyy] (red), and Im[Gzz] (blue), 
normalized to the free space value Im[Go] (dashed line). The 
number of unit cells is A/" = 81. Cesium D2-\me frequency 
va = 352 THz is centered in the shaded area. The vertical 
dotted line marks the light line, beyond which all decay chan- 
nels are lossy, (b, c) show 7tot (black curves), normalized to 
the free space value 70 (dotted line), in the frequency range 
marked by the shaded area in (a). The solid (dashed) curve is 
evaluated using 81(61) unit cells. The atomic spin is aligned 
to the X-axis, with the spin projection quantum number (b) 
mp = and (c) mp = 5. 

in jijj occur for the initial excited state 6Ps/2^F' = 
5,m^ = 0, while jid is suppressed for the initial state 
= 5. This is because the probe mode predominantly 
supports TT-polarization and hence Atuf = 0. Coupling 
between states with Atuf 7^ is small. Of course, ad- 
ditional guided modes can contribute to 7tot7 as is evi- 
denced for the rrip = 5 state due to field polarizations 
perpendicular to the atomic spin, such as z {x) for the 
single (double) beam(s) in Fig. 4(a) (5(a)). 

From 7tot(^d) and an analytic model of coupling to the 
guided- mode near the X-point, we estimate the contri- 
butions of jiD and 7' to 7tot = liD + 7' near the largest 
resonances in Figs. Qb) andjsjb) [47 . For rup = and 
A/" = 81, we find that Jid/jo — 15 and 7V70 — 1.2 for 
the single nanobeam, while Jid/jo — 21 and 7V70 — 1-0 
for the double nanobeams [52^. Here, 7o/27r = 5.2MHz, 
the free-space Cs decay rate. 

The ratios jiD/jtot and Jid/j' serve as metrics for 
the strength of atom-photon interactions for our ID pho- 
tonic crystals. For example, the resonant reflectivity tq 
of a trapped atom for the probe field should scale as 
^0 = liD/ltot El US], which for the double nanobeams 



FIG. 5: As in Fig. ^ now for the double nanobeam at r = 0. 
(a) shows the diagonal components of the Green's tensor, 
Im[Gxx\ (solid black), Im[Gyy\ (red), and Im[Gzz\ (blue), 
normalized to /m[G'o] for free space (dashed line). = 81 
unit cells. Cesium Z)2-line frequency va — 352 THz is cen- 
tered in the shaded area. The vertical dotted line marks the 
light line, beyond which all decay channels are lossy, (b, c) 
show 7tot (red curves), normalized to the free space value 70 
(dotted line), in the frequency range marked by the shaded 
area in (a). The solid (dashed) curve is evaluated using 81(61) 
unit cells. The atomic spin is aligned to the ?/-axis, with the 
spin projection quantum number (b) m'p — and (c) m'p — 5. 

leads to ro — 0.95. For a cavity QED system with 
one 'impurity' atom surrounded by Na 'mirror' atoms 
along a lZ)-lattice fT6|, the ratio of the coherent coupling 
rate gi = \^7^ji]j/2 to the effective dissipative rate 7' 
would exceed unity even for Na = 1 atom. For conven- 
tional cavity QED, we estimate a vacuum Rabi frequency 
^ 27r X 6GHz for an atom trapped in the ID-photon 
crystal waveguides studied here, making the reasonable 
assumption of a cavity formed from N ^ 10 unit cells. 

Certainly there are challenges to the implementation 
of our designs for trapping atoms in ID photonic crys- 
tal waveguides with strong single-photon interactions, in- 
cluding atom loading into the small trap volume and light 
scattering from device imperfections. We are working 
to address these issues by numerical simulation, device 
fabrication, and cold-atom experiments with nanoscopic 
structures. Our efforts are motivated by the predic- 
tion '^in/ltot ^ 0-9 in Figs. [4j [sj which is unprece- 
dented in AMO physics and which could create new scien- 
tific opportunities (e.g., quantum many-body physics for 
ID atom chains with photon-mediated interactions, and 
high-precision studies of vacuum forces). Moreover, our 
double nanobeam structure provides proof-of-principle 
for a promising new concept that combines optical and 
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vacuum forces to form stable traps for neutral atoms in 
dielectric nanostructures. 
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SUPPLEMENTAL MATERIAL 



CALCULATION OF -Hot 



CALCULATION OF CASIMIR-POLDER 
POTENTIALS 

The Casimir- Polder potential t/cp(ra) in Ref. [T is 
calculated from the following integral |2J: 



/7cp(ra) 



27r 



where Tr[.] denotes the trace, is the dynamic po- 
larizability tensor of ground-state Cesium atom, and 
G5c(ra,ra,a;) = G(ra,ra,cj) - Go(ra,ra,a;) is the scat- 
tering Green's tensor, that is, the Green's tensor G 
subtracted by the vacuum contribution Gq evaluated 
at atomic location r^; ^ttH is Planck's constant, and 
/io is the vacuum permeability. The Green's tensor 
is the solution to the Maxwell equation V x V x 

— ^e(r,a;) G(r, r',cj) = 1(5*^^^ (r — r'), corresponding to 

the electric field response to a point dipole current source. 
6(r, u) is the dielectric function, and I is the unity tensor. 

We employ finite-difference-time-domain (FDTD) cal- 
culations [3 to solve numerically for the Green's tensors 
of our structures. The integral of Eq. ([T]) is evaluated by 
adapting a procedure established in Ref. [4 and by using 
a deformed contour uj{^) = ^^^1 + ia/^ in the upper half 
of the complex frequency plane, parametrized by a real 
number ^ > and a constant a > 0. As explained in Ref. 
[4, this is equivalent to solving the Green's tensor at real 
frequencies ^ with a fictitious global conductivity applied 
to the dielectric function e'(r, ^) = (1 + zcr/^)e(r, ^). The 
integration can then be performed in the time domain 
(via the convolution theorem) and converges quickly due 
to fast decay from <j. 

Specifically, Eq. ([T]) is numerically evaluated using 
Ucp{ra) = 7^ / dtlm [g^,u{-t)] • Esc,^(ra, t), (2) 

where Esc,iy(ra,t) is the (real) electric field generated 
by a point dipole current source J = 5(t)xj, {x^ = 
x^y^z) located at the position and scattered by a 
structure with a dielectric function e'(r, ^) [5]. Here, 
the indices ji and v are repeated for summation con- 
vention, g^u{t) is the Fourier transform of g^viX) — 
-i^y/l + f (l + I) QiOalA^iO), and 6(0 is the 
Heaviside step function. For the initial calculations in 
Ref. [l] , we take the dielectric constant e to be frequency 
independent, ^(r,^^) e(r). 



To determine the total spontaneous decay rate 7tot for 
an atom in our structures, we also solve for the classical 
Green's tensors and evaluate 7tot via [H [6l [7] 



7tot 



— — ^Im 
a 



{0} 



(3) 



where Dj = ({0}|d'''|j)(j|d|{0}) is the dipole matrix ele- 
ment between the ground state manifold and the excited 
state j, and ujj is the transition frequency. The total de- 
cay rate 7tot = 7id + includes the decay rate 71D to a 
guided mode of interest as well as the rate 7' to all other 
modes of the structure, including lossy modes. As dis- 
cussed below, the contributions of 710,7' to 7tot can be 
estimated from the global frequency dependence 7tot(^)- 
To obtain jtoti^)^ we evaluate the Green's tensor for 
the real dielectric function e(r) using the FDTD method, 
followed by a discrete Fourier analysis. 

VALIDATION 

To validate our numerical procedures, we have per- 
formed calculations of Ucp for several geometries where 
analytical solutions are available, including an atom near 
an infinite dielectric or metal half-space [8 and an atom 
located above an infinite dielectric grating [9 , and found 
excellent agreement between our simulations and the ex- 
act results. 

We have validated our calculations of 7tot for the cases 
of an (atomic) dipole near an infinite dielectric, metallic 
parallel plates, a nanofiber [10 , and 2D-photonic band- 
gap microcavities [H] . 

GUIDED MODE RESONANCES 

For our structure with an infinite number of unit 
cells, a guided mode (denoted by A) contribution to 
the imaginary part of the Green's tensor can be calcu- 
lated as [7l|12|, lm[G^j^{ra,ra,uj)] = ac^ux{ra] kj,) 
u^(ra; kx)/2ujv\^ when the frequency uo intersects the fre- 
quency band uox at a propagation constant below the 
light line. Here, UA(r; kx) is the orthonormal mode func- 
tion, and vx is the group velocity, both available via nu- 
merical calculations [13l[14]. As we scan the frequency cj, 
Im [G^j^] increases monotonically and diverg es as kx ap- 
proaches the X-point, where vx 0. The guided mode 
Green's tensor vanishes when uj lies beyond the frequency 
of the band edge. 

Based on this analysis, we can evaluate the decay rate 
7^^^(co') into the designated probe mode for an infinite 
structure, and compare it with the heights of resonant 
features in 7tot('^) for finite structures with different 
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numbers of unit cells, as shown in Figs. 4(b) and 5(b) 
of Ref. [1]. Indeed, the actual 71D of a finite-size struc- 
ture must deviate from 7^^^ due to boundary conditions 
that transform a continuous spectrum into a discrete set 
of resonant peaks [15 , as shown in Figs. 4(b) and 5(b) 
[1]. When the number of unit cells is increased in our 
calculation of 71 d over the range N = 11 toA^ = 81, we 
find that the frequencies uj'^^^ of the resonant peaks shift 
in position and the peaks change height. As documented 
by the black diamonds in Figs. 1(b), 1(f) of Ref. p. , the 
cj^^^ arise from the discrete set of propagation constants 
/ci^^ Tin/aN imposed by the boundary conditions for 
the finite structures with n either even or odd. 

The peaks in 7tot at the set of frequencies uj^'^^ build 
up on top of a fairly constant background within the fre- 
quency range displayed in Figs. 4(b) and 5(b) of Ref. 
|1 . We assume that this background represents the con- 
tribution of 7' to 7tot, and subtract the background to 
estimate 71D. The resulting form for 7id(^) consists of a 
set of resonant peaks whose heights at discrete uj^'^^ qual- 
itatively map out 7^^^ calculated for the infinite struc- 
ture, with the maximum peak height for 71D occurring 
for the peak closest to the band edge. Moving further 
away from the band edge, we find that our numerical es- 
timate of 7tot(^) — asymptotes to the calculated 
value of 7^^^ {uo) reasonably well. 

Estimation of 71D and 7^ 

From the previous discussions, we identify that the 
decay rate into other modes 7' can be read off from 
the broad background in 7tot- Specifically, we estimate 
7' = 7tot(^0 ^ frequency uj' just across the band edge 
and away from any resonant peak for a guided mode. 
The decay rate into the probe mode 71D can then be 
estimated using 71D = 7tot - • 

For the single-beam structure and the atomic spin ori- 
entation shown in Fig. 4(b) [1 , we find a peak total decay 
rate 7tot/7o ^ 15 and a background level 7V70 ~ 1-2 
near the Cesium D2-\hie frequency va = 352 THz. 
We estimate the coupling to the resonant probe mode 
7iD = 7tot — 7' ~ 1470. For the double-beam struc- 
ture and spin orientation shown in Fig. 5(b), we find 
7tot/7o ~ 22, 7770 ^ 1, and, therefore, 710/70 ^ 21. 
Here, 7o/27r = 5.2MHz is the free-space (vacuum) decay 
rate for the D2 line. 

EFFECTIVE AREA AND MODE VOLUME FOR 
PROBE 

Both the single and double nanobeam structures in 
Ref. ^ lead to atom localization in a region of large am- 



plitude for the probe field. One measure of the strength 
of the atom-field coupling is the effective mode volume 
Vm per unit cell, where 

^m = j e(r)|E(r)|2dV/e(r^,n)|E(r^,n)|'. (4) 

Here the integration is carried out over the volume of a 
unit cell. That is, the integration domain along propa- 
gation direction x extends over the distance a (i.e., the 
lattice constant), while in the transverse y^z directions, 
the integration domain is from —00 to +00. 

For the single nanobeam, the probe mode has a global 
maximum at r^nin = and an effective mode volume 
Vjn ~ 0.13 jiia^ . For the double nanobeams, the probe 
mode has a saddle-like intensity distribution around 
i*min = 0, resulting in Vm ~ 0.11 jiia^ for a unit cell. 

* These authors contributed equally to this work. 
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